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Previous theoretical, along with early simulation and experimental, studies have indicated that 
particles with a short-ranged attraction exhibit a range of new dynamical arrest phenomena. These 
include very pronounced reentrance in the dynamical arrest curve, a logarithmic singularity in the 
density correlation functions, and the existence of 'attractive' and 'repulsive' glasses. Here we 
carry out extensive molecular dynamics calculations on dense systems interacting via a square-well 
potential. This is one of the simplest systems with the required properties, and may be regarded as 
canonical for interpreting the phase diagram, and now also the dynamical arrest. We confirm the 
theoretical predictions for re-entrance, logarithmic singularity, and give the first direct evidence of 
the coexistence, independent of theory, of the two coexisting glasses. We now regard the previous 
predictions of these phenomena as having been established. 



PACS numbers: 82.70Dd, 64.70Pf, 83.10Rs. 
I. INTRODUCTION 

Recently, there has emerged a series of remarkable re- 
sults involving the dynamical arrest of particles with in- 
teraction potentials that are short compared to the size of 
the repulsive core [[!]-[)) . The most important predictions 
are as follows jl(| . The curve at which the fluid phase is 
arrested is reentrant in the regime where attractions and 
repulsions compete. In the vicinity of this reentrance, 
and in the arrested regime, there exists two distinct ar- 
rested states that may, under some conditions, coexist. 
These two states differ in their long-term dynamics. This 
coexistence terminates smoothly, the dynamics of the two 
arrested states becoming identical at a particular point 
at which distinctive dynamics is expected. It is impor- 
tant to note that all these phenomena are predicted to 
occur irrespective of the shape of the potential (explicit 
results exist for square well or hard-core Yukawa interac- 
tion models ) ■ The essential feature is only that the 
range of the attraction in comparison to the hard core be 
extremely small. 

The results described above were first deduced from 
the mode-coupling-theory (MCT) [0,0, and we may 
note that the conclusions do not depend on de- 
tails of the approximation for input static structure, 
Percus-Yevick, Mean Spherical Approximation and Self- 
Consistent Ornstein-Zernike Approximation, all leading 
to essentially the same picture. Within MCT the merg- 
ing of the two arrested states is predicted to be a higher 
order dynamical singularity, and density correlations in 
its vicinity are expected to obey a highly distinctive log- 



arithmic relaxation j|,|l3| , in contrast to the conventional 
ergodic-non ergodic transition where a two-step process 
takes place @|4|. 

More recently this decay has been observed also in sev- 
eral experimental systems p^| , p^ ,p| , thus giving the first 
evidence that these singularities do exist in nature. 

Two recent numerical works [p||l7j] have focused on 
the new dynamical features characterizing attractive col- 
loidal systems. In Ref. ||, colloidal interactions were 
modeled with a potential chosen in such a way to avoid 
undesired effects such as liquid-gas separation at low den- 
sities. Polydispersity was included to prevent crystalliza- 
tion at high densities. In Ref. Jl7j Molecular Dynamics 
simulations of a square well (SW) system with very nar- 
row range of attraction have been performed. In both 
cases, the results have been shown to be in excellent 
agreement with MCT predictions. The focus on the sim- 
ple square well potential, makes direct contact both with 
the available theoretical results and with experimental 
systems, being the interactions completely controllable 
and not dependent on external parameters. For the SW 
model, theoretical solution is also available. Indeed, the 
equilibrium properties of square- well potential have been 
studied for many years, and in many ways is regarded as 
the simplest physical model that exhibits all of the essen- 
tial features of the short-ranged systems |l^-|2^| . Besides 
this, it is the model for which the theory of the dynami- 
cal arrest transition described above has been developed 
and for which, therefore, detailed comparison may 
be made between theory and simulation. It may therefore 
be regarded as canonical in this arena of study. 

In Ref. [p~7fl , we examined a one-component square well 
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system. We found that, no matter the sharp intervening 
of crystallization which effectively prevented the system 
to approach very closely the glass transition, it was pos- 
sible to have a clear picture of the reentrant shape of the 
glass curve in the temperature-density plane. This was 
achieved by plotting iso-diffusivity curves, and examin- 
ing their trends when approaching the limit D — ► 0. The 
shape of the liquid-glass line found was in good agree- 
ment with the previous MCT 0] predictions. 

In this paper, we report an extensive numerical study 
of a two-component binary mixture with interactions 
modeled by a very narrow SW potential. This mixture 
appears to be a logical extension of the one-component 
system for which crystallization is effectively avoided, 
by means of the geometrical rearrangements allowed by 
asymmetry in diameters of different particle species. The 
choice of a two-component system makes possible to 
extend the range of iso-diffusivity curves of almost 4 
decades, as well as the range of studied packing fractions 
from approximately 0.57 for one-component to 0.62 for 
binary systems. One notes that at that highest density, 
the reentrance is so pronounced that an equivalent hard 
sphere system would be sufficiently dense to approach 
the 'random-close-packed' limit, though the definition of 
that concept has its own limitations p3]. 

The extension to a binary system does not only make 
more evident the reentrant behaviour of the glass transi- 
tion curve, but also allows a deeper study of the dynam- 
ics of the system. By this we mean that it turns out to 
be possible to reach very closely the ideal glass line, in- 
tended as the D — > locus of the (cf>, T)-plane, and thus 
to feel quite clearly the presence of higher order MCT 
singularities. We know from the theoretical work that 
the 3%-case of the square well model possesses indeed an 
^-singularity ^,||, corresponding to the end-point of 
the repulsive-attractive glass-glass transition. Though, 
of course, this singularity lies inside the glassy region, its 
presence is signaled by the characteristic logarithmic de- 
cay of density correlators mentioned above, already in the 
liquid region, when going sufficiently close to the neigh- 
bouring glass boundary. 

We report in this paper a characteristic behaviour of 
the density correlators, near the ideal glass transition, 
which combines features of the typical Ai singularity, i.e. 
simply ergodic to non-ergodic transition with two-step 
power-law relaxation, and of higher order singularities 
associated with logarithmic decay. This produces, in a 
certain region of the (<fi, T)-planc, a new behaviour which 
cannot be described by the asymptotic predictions for ei- 
ther of the two behaviours, deriving by a competition of 
the two types of MCT solutions. The same kind of in- 
teresting result is found in the mean square displacement 
(MSD). This subtle interplay between different singular- 
ities is present also in theory, and its manifestation in 
simulation may be regarded as support for the theory, 
rather than an inconvenience. This will allow us to iden- 



tify and localize quite clearly a genuine MCT higher order 
singularity in a realistic model. 

II. SIMULATION AND THEORY 

We study binary mixtures of SW spheres. In partic- 
ular, we focus on samples of a 50% — 50% mixture of 
N = 700 spheres in a cubic box, with diameter ratio be- 
tween the two species equal to 1.2. Thus, the smaller 
particles (B-type) diameter is <jb — 1, and both A and 
B particles have unit mass. Both species interact with a 
square well potential with ratio between potential range 
and particle diameter equal to 3%. This corresponds to 
one of the cases studied theoretically within MCT, that 
clearly possesses all the main phenomena |,|,||. The 
3%-ratio has been chosen also for interactions between 
particles of different species, i.e. we consider 

Vij (r) = 00 r VJ < (Tij 

Vij(r) = -V.Q (Tij < Tij < (Tij + Ay 

Vij(r) = nj > (Tij +A i:j (1) 

with eij = Ay*/(er.y + Aij) = 0.03, i,j — A,B and we 
use the conventional notation for which, for example, 
&aa = &A and gab — (&A + o\b)/2. Temperature T 
is measured in units of energy, i.e. fcs = 1 and thus, 
for example, T = 1 corresponds to the system having a 
thermal energy equal to the well-depth, while the pack- 
ing fraction is defined as <j> = (pa&a + Pb^b) ' n /®i where 
Pi = Ni/L 3 , L being the box size and iVj the number of 
particles for each species. 

Initial configurations for each density were chosen at 
random. Particles were separated in successive steps, 
with more particular care the higher the density of in- 
terest, to implement the hard core repulsion. When sep- 
aration was ensured, the attractive well was added. To 
reach the temperature of study, the configuration so pre- 
pared was then left to evolve with a thermostat of con- 
stant thermal coefficient for a period of time sufficient to 
equilibrate at that temperature. We estimate the equi- 
libration time as the time at which the density correla- 
tion function of the slowest collective mode (i.e. at the 
structure factor peak) has decayed to zero. After this 
equilibration, the configuration was left to run at con- 
stant energy for a time dependent on the slowness of the 
dynamics, for a time covering at least 10 equilibration 
times. 

Simulation time for each species is measured in units 
of <7j • (m/ito) 1 / 2 . Standard MD algorithm has been im- 
plemented for particles interacting with SW potentials 
p6| . Between collisions, particles move along straight 
lines with constant velocities. When the distance be- 
tween the particles becomes equal to the distance for 
which V(r) has a discontinuity, the velocities of the inter- 
acting particles instantaneously change. The algorithm 
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calculates the shortest collision time in the system and 
propagate the trajectory from one collision to the next 
one. Calculations of the next collision time are optimized 
by dividing the system in small subsystems, so that col- 
lision time are computed only between particles in the 
neighboring subsystems. 

We studied eight isothermal cuts of the phase dia- 
grams, with temperature varying between 2.0 and 0.3 in 
the large packing fraction region, i.e. <j) > 0.5. In addi- 
tion, we examined the hard spheres case, where no attrac- 
tive interactions are present. For each considered con- 
figuration, we first studied the thermal history to check 
that, effectively, it maintains itself at the required tem- 
perature within fluctuations and that the total energy 
remains constant. 

Of each studied configuration, we considered the time- 
dependent density correlation functions for different q- 
vectors to make direct comparison with the behaviour 
predicted by MCT. The correlators are defined as, 



m-fa~hi 1] (t/to)- a + hW(t/t )- 2a 
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Sij (q, t) = (Pi(-q,t)pj(q,0)) 
Sij(q) (pi(-q,0)pj(q,0)) 



(2) 



where p(q, t) = exp (iq • r;(t)) are the density vari- 
ables for each species and Sij(q) are the partial static 
structure factors of the system. The <fiij(q,t) are the 
fundamental quantities of interest in the Mode Coupling 
theory, which consists of writing a set of generalised 
Langevin equations, that can be closed within certain 
approximations [0,0,[27|,|8i . 



The correlators <f>ij(q,t) have been calculated averag- 
ing over several independent configurations and over up 
to 100 different wave-vectors with the same modulus, to 
obtain a good statistical sample. Of course, these quan- 
tities show an interesting behaviour where the dynamics 
of the system gets slower, i.e. in the supercooled regime. 
Indeed, the relaxation process starts to show a separation 
in two time-scales, which originates the typical two-step 
relaxation scenario near the (ideal) glass transition. A 
first relaxation process, the so called /3-relaxation, oc- 
curs at short times, due to particles exploring the cage 
formed by their neighbours, while a second process, the 
a-relaxation, that accounts for the restoration of the er- 
godicity due to breaking of the cages, occurs at larger 
and larger time scales, due to the slowing down of the dy- 
namics, forming a longer and longer plateau region. At 
the ideal glass transition, the time of the a-relaxation di- 
verges, and the correlators do not relax anymore, thus re- 
maining at this plateau value. This is defined as the non- 
ergodicity parameter fij(q) = lim^oo <pij(q, t), which 
jumps discontinuously from zero to a finite (critical) value 
fij(q), signaling the occurrence of an ergodic (fluid) to 
non-ergodic (glass) transition. 

The two-step relaxation is well described by MCT, 
through an asymptotic study of the correlators near the 
ideal glass solutions. The approach to the plateau is de- 
scribed by a power law, regulated by the exponent a, i.e. 



with to the microscopic time, while the departure from 
the plateau, i.e. the start of the a-process, is expressed 
in terms of another power law, regulated by the exponent 
b, 



2 b 



(4) 



with t the characteristic time of the relaxation. The 
exponents a and b are related to each other with an al- 
gebraic relation, and are independent of the particular 
q-vector considered, h q x ^ and h q are called critical am- 
plitudes (T||. 

On the other hand, the a— relaxation process can be 
also well described by a stretched exponential, i.e. 



dp q {t) = A qC *v[-(t/T q y«] 



(•5) 



where the amplitude A q determines the plateau value, 
and the exponent (3 q is always less than 1. 

We can thus fit the correlators, for each (/-value con- 
sidered, both in terms of the MCT prediction and of the 
stretched exponential, and find an estimate of the non- 
ergodicity factor f q as a function of q. The shape of this 
quantity can be indicative of the formation of different 
types of glasses, either attractive or repulsive dominated 

However, what we have described so far is typical for 
Ai singularities. These correspond to the simplest non- 
trivial solutions for the non-ergodicity parameter MCT 
equations, and for example in the hard sphere model only 
this type of singularity can arise. This is due to the 
fact that the only control parameter of the model is the 
packing fraction. When the number of control parame- 
ters increases, higher order singularities may occur. For 
a square well model it was shown that singularities of 
types A3 and A4 are present within the theory Q when 
the width of the well becomes much smaller than the 
hard core radius. In the proximity of such singularities, 
the asymptotic behaviour for the density correlators is 
different from the one we have seen so far 0,[l3| . In par- 
ticular, we have for the leading contribution a logarithmic 
behaviour, 
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Another main focus of our study was to evaluate the 
MSD (r 2 (t)) of particles with respect to their initial posi- 
tions. We considered a configuration to have sufficiently 
evolved when particles on average have traveled further 
than a few diameters. 

Typically, the behaviour of the MSD at short times, 
follows the simple law (r 2 (t)) ~ t 2 , which accounts for 
the ballistic motion of the particles, i.e. particles move 
freely without collisions. At later times, particles start to 
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feel the presence of each other and there is crossover to 
the diffusive regime, i.e. (r 2 (t)) ~ t. The proportional- 
ity constant of this relation defines the diffusivity of the 
system, via the celebrated Einstein relation pjJ, 



lim 

t — >oc 



(r\t)) 



t 



6D 



(J) 



Thus, from evaluating the long time limit of the MSD, 
we determine the diffusion coefficient D for each state 
point. 

Though this general behaviour is preserved, when the 
dynamics becomes slower, a transient region between 
short and late times emerges. Of course, the duration 
of this transient region increases the slower the dynam- 
ics. This phenomenon is the correspondent behaviour for 
the MSD of the separation of the two time scales that we 
have seen in the correlation functions. It also reflects the 
formation of cages in which particles get trapped, so that 
diffusion becomes more difficult. The crossover region 
consists generally in the development of a plateau also for 
the MSD. At the ideal glass transition, the a— relaxation 
time would diverge, and diffusion from the plateau would 
not occur even at infinite times, in complete analogy with 
the correlators behaviour. Thus, the height of the plateau 
represents the localization length of particles in the ar- 
rested state, i.e. the size of the cages of the glass. Of 
course, in simulations, only finite times can be explored 
and the position of the ideal glass transition can be only 
extrapolated by data. 

As for the correlators, the presence of higher order sin- 
gularities may affect the plateau region of MSD, giving 
rise to peculiar new behaviour. 

MCT predicts a power law decrease for the diffusiv- 
ity on approaching the ideal glass transition. Along an 
isotherm 
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(8) 



where 4>c is the value of the packing fraction at the transi- 
tion ('critical' value), i.e. the value where the diffusivity 
would drop to zero for the considered temperature. The 
exponent 7 is completely determined by the theory in 
terms of the exponents a and 6, via the simple relation 
7 = l/(2a) + 1/(26) @. It is also related to the so-called 
exponent parameter A by an analytical relation. This pa- 
rameter is crucial in determining the presence of higher 
order singularities 0, in particular it tends to the value 
A = 1 at an A3-point, while for a simple Ai singularity 
is always less than 1. 

So, in principle, from fitting the diffusivity behaviour 
at constant temperature with Eq. |], one can determine 
the exponent 7, and from this, also a and b. From a 
fitting procedure based on Eq. ^, the non-ergodicity pa- 
rameter can be calculated consistently. However, close 
to higher order singularities, we do not expect that this 
behaviour is generally preserved, as the logarithmic be- 
haviour in the correlators and the transient region in the 



MSD intervene. Indeed, in these conditions, the expo- 
nent b tends to zero (thus originating the logarithmic 
behaviour) and, from the relation between b and 7, one 
can see that 7 would go to 00, and Eq. |^ cannot then de- 
scribe the arrest. Thus, the region of validity of asymp- 
totic predictions may shrink significantly close to higher 
order singularities. 



III. RESULTS: THE OVERALL PICTURE 

The considered mixture represents a natural extension 
of the mono-disperse system studied so far jl7|]. Indeed, 
the small asymmetry in diameter does not produce sig- 
nificant change in the dynamics of the two cases, and, on 
the contrary, allows to reach much larger packing frac- 
tions with no sign of crystalline order. The first of the 
two sentences can be explained by looking at Figure [I]. 
Here, we compare the behaviour of density correlators at 
a correspondent point on the control parameter space of 
the system. For the mono-disperse case we are almost 
at the most reentrant point before crystallization takes 
place, i.e. <j> = 0-57 and T — 0.75. We can clearly see 
that, nonetheless, dynamics does not appear to be par- 
ticularly slowed down. To make the comparison with the 
binary case, we are considering the total density correla- 
tion function for the species 1 at the g-vector correspond- 
ing approximately to the first peak of the static structure 
factor, and quantities have been rescaled in order to com- 
pare particles with equal diameters. It is evident that the 
dynamical behaviours are very close, thus in this sense, 
we can think of doing an extension of the one-component 
work. 

In the following, we will focus on properties of particles 
of type A. Thus all quantities reported without label, will 
refer to them. This choice derives from the fact that we 
do not expect substantial differences in the behaviour of 
the two species, due to the small amount of asymmetry 
in their sizes. 

We start by comparing the T and cf> dependence of the 
diffusion coefficient. The diffusion coefficients can be nor- 
malised with respect to the factor D % = ai y^T/m, which 
takes into account the T-dependence of the microscopic 
time. This ensures that the difference in the average ve- 
locities due to the temperature is eliminated, and the 
diffusion can be considered to be comparable between 
different temperatures. 

A plot of the Da/ Dq in function of packing fraction, 
along the considered isotherms, is shown in Fig. || In 
the present work, we focus our attention mainly on the 
high densities regime, i.e. <f> > 0.5. The behaviour of 
the diffusivities presents many similarities with the case 
of a mono-disperse sample of SW spheres [^|. However, 
striking novel features appear due to the fact that, for the 
chosen binary system, it is possible to reach diffusivities 
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4 orders of magnitude smaller than in the mono-disperse 
case, as well as much larger packing fractions with no 
sign of crystallization. 

We present results of normalised diffusivities varying 
roughly between 10 -2 and 10~ 6 , while the mono-disperse 
system could only reach values up to slightly above 10~ 2 , 
due to the intervening of crystallization |l7| . We remark 
that the lowest diffusivity values were imposed by compu- 
tational times, and not as in the mono-disperse by crys- 
tallization. Also, the attractive binary system is able to 
occupy effectively a larger amount of available volume, 
thus reaching liquid states up to a packing of about 0.62. 
On the other hand, the hard spheres case reaches compa- 
rable values of diffusivities at a packing fraction of about 
4> = 0.585 (see Fig. f|), a value close to the one experi- 
mentally established for the one-component hard sphere 
case §§. 

Examining the figure, it is evident that the behaviour 
of the diffusivity is driven by two competing mechanisms. 
Upon decreasing temperature starting from the highest 
value, the presence of the repulsive core, initially domi- 
nant, enters in competition with the attractive interac- 
tions. This is manifested in the diffusion getting larger, 
at the correspondent packing fraction, by decreasing tem- 
perature. In other words, the system reaches the same 
diffusivity at larger and larger packing fractions. This 
is due mainly to geometrical rearrangement of particles, 
as the temperature is lowered, i.e. particles tend to get 
closer and, consequently, more free space for diffusion 
gets available to the system. However, when the temper- 
ature becomes small enough, i.e. effectively lower than 
the energy scale of the square well, attractions become 
dominant and thus, diffusion becomes slower again be- 
cause particles tend to remain within each other shell of 
attraction. 

We know from theoretical calculations within MCT 
that this phenomenon is typical of very narrow attrac- 
tive potentials, both for square well interactions jl and 
for hard-core Yukawa ||. Indeed, the 3%-choice for 
the range of the attractive well in our simulation cor- 
responds to ensure that competition between attrac- 
tion and repulsion is particularly sharp. This can be 
explained in terms of cages, i.e. when the attractive 
range is not small enough, which means few percents 
of diameter, there is not much difference between cages 
formed by neighbouring particles at high or low temper- 
atures. On the other hand, a very localised attraction 
can effectively change the shape of the cages, by stick- 
ing particles within the well-distance A. This produces 
the larger diffusions observed at intermediate temper- 
atures, when the two mechanisms almost balance each 
other. Similarly, extremely short-ranged attractions pro- 
duce a solid-solid iso-structural transition, between an 
attractive-dominated and a repulsive-dominated crystal 
pl[ . This has been correlated with the glass-glass tran- 
sition predicted by MCT in a recent work 0). Of course, 



it would be interesting to extend the simulations to dif- 
ferent values of the range of the potential to confirm the 
width dependence of the anomalous behaviour. 

Observing more carefully Figure ||, the two mecha- 
nisms of diffusion produce two different trends of be- 
haviour for the plotted curves. Indeed, for T > 0.6 
the curves present a quite dramatic decrease of diffusiv- 
ity, while for smaller temperatures the same decrease of 
about 4 orders of magnitude occurs on a much wider 
range of packing fractions (for example comparing T — 
2.0 and T = 0.4 such range almost doubles). 

We also note that, if we plot the bare diffusion coef- 
ficients, as evaluated from the fit of the MSD, without 
normalising by D°, the reentrant behaviour is preserved. 
This can be an advantage, from an experimental point 
of view, because it would allow to observe at the same 
packing fraction for various temperatures the diffusivity 
firstly increasing then decreasing again, without having 
to include the thermal factor. We plot in Fig. || the not- 
normalised diffusivities at fixed packing fractions, vary- 
ing the temperature. The appearance of a maximum, 
sharper with increasing density, is indeed another mani- 
festation of the reentrance. Extracting the values of max- 
imum diffusivity, it is possible to draw a 'maximum diffu- 
sion' line in the phase diagram. This could be interesting 
as, for example, in a complex system like water, this line 
appears to play a very important role in the understand- 
ing of the 'metastable' part of the water phase diagram 
& 

To make contact with the theoretical results for the 
ideal glass transition, we have extrapolated curves of nor- 
malised iso-diffusivities, as for the mono-disperse case 
p7[ and represented them in Figure ||. Of course, the 
limit D — > would correspond to the ideal glass line, as 
calculated by MCT. We report the curves for normalised 
values varying between 5 • 10 -3 and 5 ■ 10 -6 , and in the 
inset we present for comparison the ideal glass line as 
predicted by MCT for a SW one-component system. 

It is interesting to comment on the behaviour of the iso- 
diffusivity curves. Depending on which diffusivity value 
is chosen, the most reentrant point, i.e. the (</>, T)-point 
characterized by largest packing fraction with that diffu- 
sivity, changes. In the considered range of diffusivities, 
its temperature varies from 0.75, which also corresponds 
to the most reentrant point for the mono-disperse case, 
to 0.5. However, data of Fig. || allow to say that, in the 
limit D — > 0, such point will be found at a finite tem- 
perature between 0.4 and 0.5, since the reached packing 
fractions are so large that it would not be possible, within 
the trend, to go much beyond. 

The reentrant behaviour, present in Mode Coupling 
calculations, is then confirmed by simulation. It is clear 
that, since MCT underestimates the effect of packing, we 
should not expect, a perfect quantitative agreement be- 
tween theory and numerical results. Indeed, for example 
for the simple hard sphere case, the MCT critical glass 
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transition packing fraction is 4> ~ 0.516 |33| whereas in 
experiments this has been shown to be (j> ~ 0.58 f30"| . 
On the other hand, the same experiments have shown 
that more accurate predictions can be expected for the 
behaviour of dynamical quantities such as the density- 
density correlation functions and the MSD. 

IV. RESULTS: DYNAMICS ALONG ISOTHERMS 

We start by examining the results for the isotherm 
T = 2.0. Here particles have sufficient thermal energy 
to escape the attractions and the resulting dynamics 
appears to be quite similar to that expected for ordi- 
nary hard spheres. However, the effect of the attraction, 
thought not changing the general behaviour, is to en- 
large the liquid part of the phase diagram toward pack- 
ing fractions already quite larger than the typical (one- 
component) hard sphere value, i.e. 0.58. At this temper- 
ature, indeed, the system behaves as a fluid at least up 
to a packing fraction <f> — 0.595, where the time limit of 
our simulation is reached. Close to this limit value, the 
diffusivity decreases almost two orders of magnitude for 
a variation of 1% in packing fraction. 

In Figure |s](a) we report the time-dependent density- 
density correlators at increasing packing fraction, up to 
the closest to the ideal glass transition. The dramatic 
decrease in diffusivity is reflected in the behaviour of the 
correlators by the formation of the typical two step relax- 
ation process near the arrest, described above. Similar 
behaviour can be also observed in the behaviour of the 
MSD displayed in Fig. ^(b) in unity of a 2 . We note, 
however, that the height of the plateau of approximately 
10% of particles diameter is consistent with Lindemann's 
melting criterion |34J . 

The diffusivity data can be fitted with the MCT power- 
law behaviour (Eq. ||) . This holds sufficiently close to the 
ideal glass transition, thus we considered relevant for the 
fit only those points for which a clear a-relaxation process 
was evident. Doing so, we found 7 ~ 3. 6 ±0.4. This value 
is already much larger than the typical one-component 
hard sphere value predicted by MCT, i.e. j HS = 2.58 
p3| . In the case of our particular binary mixture, in 
the hard sphere case, whose diffusivity behaviour is also 
reported in Fig. 0, we found a value of 7 = 2.9 ± 0.2. 
The uncertainty on the exponent is due to the variation 
it gets when considering only the points closest to the 
transition. We note that in reference the case re- 
ported brings a value of 7 = 3.03. Thus, we are in a sit- 
uation clearly closer to a higher order singularity, being 
the correspondent other MCT exponents for this value, 
respectively b = 0.35 and A = 0.874 ± 0.03. It is perhaps 
important to stress at this point that the value of 7 ob- 
tained with this procedure can be slightly wrong due to 
the difficulty to get close enough to the ideal glass tran- 
sition with numerical simulations iBal. We will use the 



so-calculated value of b in the next paragraph to fit the 
behaviour of density correlators along the iso-diffusivity 
curve D/Dq = 5 • 10~ 6 , to give an idea of the behaviour 
of quantities of interest along the ideal glass transition 
line. 

The case T = 1.5 shows a behaviour completely anal- 
ogous to T = 2.0. We note however that the fit of the 
diffusivity with Eq. [8] gives in this case the exponents 
7 = 3.8 ± 0.5, b = 0.325 and A = 0.888. The increase in 
the value of A is expected since we are getting closer to 
the reentrant region, and consequently to the singular- 
ity. The same trend is observed for T = 1.0 with 7 still 
increasing up to about 3.96 ± 0.12, and A reaching the 
value of 0.896. 

At T = 0.75 a new interesting behaviour appears. In- 
deed, in the density correlators a logarithmic decay starts 
to emerge. As shown in Fig. |^, some state points display 
a logarithmic decay (see curve for = 0.58 in the figure) 
for almost the whole relaxation process, i.e. after the mi- 
croscopic relaxation up to complete decay. The shape of 
this logarithmic behaviour appears quite different from 
the one found in reference On the other hand, it re- 
minds quite closely the shape of MCT correlators near 
the A3 singularity for the 3% square well potential re- 
ported in Figure 11 of Ref. Q, which for comparison is 
reported in the inset in Fig. |(| Upon increasing den- 
sity, and so getting closer to the glass transition, the 
relaxation changes to the usual two-step form, clearly 
indicating a similar situation of our isothermal path to 
that indicated in the inset of Fig. 11 cited above. Thus, 
the higher order singularity dominates the dynamics at 
smaller packing fractions, but when one gets sufficiently 
close to the glass transition, a conventional A2 singularity 
is met, and this causes the restoration of the typical a- 
relaxation. By considering only those packing fractions, 
when at least the beginning of the a-relaxation can be 
observed (i.e. <j) > 0-6), we can fit the diffusivity again 
with Eq. |j, obtaining the extremely high value 7 ~ 5.1, 
corresponding to the value of b ~ 0.25 and A = 0.937. We 
cannot here estimate the error, due to the small number 
of points available (only five for three parameters for the 
fit), however, even if not so precise, it clearly indicates 
we are approaching a higher order singularity. As we 
shall see by looking at the shape of the non-ergodicity 
factors in the next paragraph, we can anticipate that we 
are meeting the Ai line along its repulsive branch, again 
as in the inset of Fig. 11 of ref. Q, but probably at a 
slightly higher temperature. 

The last observation derives from the analysis of the 
T = 0.6 isotherm. The correlators, reported in Fig- 
ure ^|(a), show an even closer behaviour to the one pre- 
dicted by MCT in @] (again Fig. 11, i.e. inset of Fig. |). 
However, at this temperature, the two competing sin- 
gularities must be so close to each other that a clear 
a-relaxation does not take place within the reach of our 
simulation, i.e. the logarithmic behaviour remains always 
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very important, and even at higher packing fractions it is 
clearly observable before the a-process takes over. The 
interesting feature emerging is that, in all the cases con- 
sidered the logarithmic behaviour does never extend for 
much more than 3 and a half decades in time. This arises 
because, in the present topology of the phase diagram, 
one is either to close to the A2 -singularities to observe 
a pure logarithm, or is too far from the glass transition 
and thus the relaxation time is generally not too large. 
Indeed, this behaviour is strongly supported again by the 
theoretical calculations in 0). 

It is to note that at this temperature we are not able 
to convincingly fit the power-law density dependence of 
the diffusivity. Indeed, if one forces the fit on the points, 
one find exponents strongly dependent on the selected 4> 
range. A possible explanation for this data sensitivity to 
4> can be found in the competition between type A2 and 
A3 or A4 dynamics. In such condition, only a compari- 
son with a full MCT solution (as opposed to an asymp- 
totic prediction) may help in rationalizing the density de- 
pendence of diffusivity. In agreement with the previous 
observations, also the MSD, represented in Figure ^(b) 
starts to show deviations from usual type B behavior. 
Indeed, a clear flat region does not appear, though dy- 
namics are significantly slow. What can be observed is 
a slight deviation from the flat region found at higher 
temperatures, which will become more and more evident 
at lower temperatures. No clear localization length can 
be found. Attraction at this temperature has become 
quite relevant. It is again a sign of very strong compe- 
tition between different singularities, between attractive 
and repulsive cages. 

Upon further decreasing temperature, we enter the 
most delicate region of the phase diagram. Indeed at 
T = 0.5, as for T — 0.6, it is not possible to find any 
MCT exponents, and the situation gets even worse in 
interpreting the behaviour of the correlators. These are 
plotted in Figure ||(a). Indeed, no clear behaviour, at all 
densities, either logarithmic or of type a, can be indi- 
viduated, and no evident plateau value ever arises. We 
note that the long-time decay of density correlators can 
be represented by a stretched exponential, but with very 
low exponents f3 q , as shown in Figure ^. In the MSD, 
reported in Figure ||(b) , the phenomenon present at the 
previous temperature, becomes more accentuated. Even 
the slowest studied state point is far from being asymp- 
totic, and the MSD presents a clear transient region. 

The case where the anomalous dynamics and the in- 
terplay between different singularities is fully displayed 
is offered by the T — 0.4 isotherm. The correlation 
functions, shown in Figure [To|(a), are rather peculiar. 
Even the long time limit is far from being rationalized 
in term of stretched exponential decay. The MSD be- 
haviour, shown in Figure ^(b) is also quite intriguing. 
The MSD transient behaviour is now evidently of a pe- 
culiar type. Indeed, for about 4 decades in time, it shows 



a dependence which can be quite accurately described by 
a power-law behaviour, i.e. 

{r 2 )~t* (9) 

We estimated via a fit x ~ 0.44. A similar behaviour has 
been found in the MCT study of polymeric systems |36) , 
for displacements varying from the typical localization 
length of hard-sphere-like cages to end-to-end distance. 
The analogy with the polymeric systems, where perma- 
nent bonds are present (in a sense close to the attractive 
cages at this very low temperature), can be a guide to a 
deeper understanding of this regime. 

The strong effects that we find at this value of temper- 
ature seem to suggest that along this isotherm the system 
approaches the closest point to the singularity, even if we 
do not know yet on which side (attractive or repulsive) 
of the glass line it will be located. To understand the 
nature of the dynamics which takes place here further 
investigations are needed, and maybe a more complete 
analysis of the correlators, and a comparison with full 
solutions of the MCT equations. 

Finally we analyze the last isotherm, corresponding at 
T = 0.3. This being a very low value for the system to 
equilibrate, data are not so clean as for the other cases, 
also because here one needs to study slower points with 
respect to the other temperatures in terms of bare dif- 
fusivities, to reach the same values of normalised ones. 
However, despite these technical difficulties, we find more 
transparent results in terms of conventional MCT inter- 
pretations, i.e. we can identify the development of a two- 
step process typical of an A2 singularity, both for the 
correlators and for the MSD than in the previous case. 

Indeed, observing the correlators in Figure [ll](a), it is 
clear that, close enough to the transition, they present 
the development of a plateau, and thus, an a-relaxation 
process, as shown in the inset of the figure for various q- 
vectors. However, this plateau is extremely high, and 
clearly indicates that we are now undoubtedly in the 
condition of attractive glasses. It is of deep interest to 
note the analogy of this behaviour with the one that has 
been found in the study of 'strong' gels |$7|, even if this 
should be inspected at lower densities also. Despite this 
clear behaviour, even at this temperature, it is not possi- 
ble to evaluate unambiguous power-law exponents from 
the diffusivity. Figure [T^(b) , represents the evolution of 
MSD at this temperature. Here, a signature of a local- 
ization length, smaller than the one in the high temper- 
ature cases, starts to develop. Indeed, an indication of 
a plateau is observable around (r 2 ) = 0.0006 • a 2 . The 
corresponding localization length is of the order of A, 
supporting the interpretation that at this temperature 
the relevant localization length has become the attractive 
well. In this respect, one can interpret the sub-diffusive 
behavior of the MSD discussed at T = 0.4 as a cross-over 
effect between the hard-sphere and attractive well differ- 
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ent localization lengths, in a similar fashion to what has 
been obtained for polymeric systems |36| . 

V. RESULTS: ALONG THE ISO-D/D CURVE 

We now focus on studying the behaviour of corre- 
lators and MSD, and other quantities along the iso- 
(normalised)-diffusivity i.e. \so-D/D , curve D/D = 
5 • 1CP 6 , shown in Fig. [|, which represents our closest 
available representation of the ideal glass transition line. 
The aim of this study is to give clear evidence of the 
existence of two distinct glassy states, attractive and re- 
pulsive. Also, it aims to connect even more closely this 
simulation to the MCT calculations, which also were per- 
formed in a similar fashion, along the ideal glass lines, in 
ref. §. 

We start by representing the behaviour of density cor- 
relators along the iso-D/Dq line in Figure [l^. The curves 
here represented, having equal diffusivity, also have the 
same normalized relaxation time. Thus, we can clearly 
see the change in the decay which takes place, upon 
decreasing temperature, from a marked a-relaxation at 
higher temperature to the extremely slow decay of T — 
0.4, passing through the intermediate regimes between 
T = 0.75 and T = 0.5. Of course, here, no evident loga- 
rithmic behaviour can be observed, due to the proximity 
to the A2 transition, as discussed in the previous para- 
graph. 

Next, we report the MSD behaviour along the line in 
Figure [l3|. Despite the larger statistical error at T = 0.3, 
we display them as an important part of the whole pic- 
ture. Thus, here, we can clearly observe the change in 
the diffusion process. The first evident thing to note is 
the big gap, of almost 2 orders of magnitude, in the lo- 
calization length of particles, which, as discussed above, 
characterises the size of cages around particles. Clearly 
this fact can be used as a justification to speak of 'attrac- 
tive' cages, opposed to normal cages, intended as simple 
occupation of the available space. In the attractive cages, 
the average distance between particles is much smaller, 
the lower the temperature. This is the first clear signal 
of different structure in the glass formation. 

Also, we can also examine more carefully the modifi- 
cation of the plateau present at higher temperatures, in 
the transient regime. Indeed, increasing the attraction, 
this tends to bend downwards until a sort of 'saturation' 
between the two competing mechanisms (attraction and 
packing) take place, corresponding to the sub-diffusive 
behaviour of T — 0.4. After this point, attractions be- 
come dominant, and the curve starts to bend upwards. 
This might suggest that, going at even lower diffusivi- 
ties, the MSD would display a similar plateau as for high 
temperatures at roughly (r 2 ) ~ 0.0007 ■ a 1 p8| , which 
means roughly a localization length of 2.6% of the par- 
ticle diameter, i.e. comparable with the width of the 



attractive well of the model, confirming our conjectures 
on the formation of attractive cages, or, to use a another 
expression, bonds. However, to gain further evidence on 
how these mechanisms really happen and evolve in the 
system, a specific study of configurations in terms of av- 
erage distance, sizes of clusters, and heterogeneities in 
general should be performed, and this is beyond the aim 
of the present work. We note that a similar figure, show- 
ing the behaviour of MSD with attraction, has appeared 
in II , but not all of these considerations could be made 
there, due to the distance from the transition. 

We now turn to evaluate the non-ergodicity factor f q 
along the iso-D/Z?o curve. To do this, we have fitted the 
density correlators at various g-vectors, and extracted 
the relevant parameters. Where possible, i.e. where 
the power-law behaviour for the diffusivity in (^|) was 
found to be valid, we used the power-law described in 
terms of the b exponent for the a-relaxation of equa- 
tion (^). Thus, for T = 2.0, we implemented the fit 
with b = 0.35, while respectively for T = 1.5, 1.0, 0.75 we 
used b — 0.32,0.31,0.25. These values have been found 
very good for the fits, always finding a \ 2 of the order 
of 10~ 4 or less. For lower temperatures, it was not pos- 
sible to use this strategy and, consequently, we had to 
use the approach of the stretched exponential in Eq ||, 
and use its amplitude as an estimate for f q . The param- 
eters of the fits, i.e. the exponent of the stretching (3 q 
and the relaxation time r g , for the considered tempera- 
tures T — 0.6, 0.5, 0.3 are reported in Fig. || to display 
their q— behaviour. Even though the stretched exponen- 
tial law is not analytically justified, it is quite established 
in the literature to use it as a fitting law for extracting 
the non-ergodicity parameter [Q. In the case of T = 0.4 
also this strategy did not work, as already discussed |39) . 

In Figure [l4| we show the so-calculated f q . Amaz- 
ingly, from T = 2.0 down to T = 0.6, they all collapse 
onto the same curves, giving a strong evidence of MCT 
predictions for the repulsive glass ||. Thus, the repul- 
sive glass is independent of temperature, and also this 
shows how the passage to attractive glass intervenes quite 
sharply. For lower temperatures, the glass becomes then 
attractive, and the non-ergodicity parameter starts to be 
modified with temperature, becoming finite also at much 
larger q- vectors [pi . Despite some errors generated by the 
stretched exponential fits at these low temperatures or 
the data noise at T — 0.3 a significant change in the shape 
and width of f q is seen between T = 0.5 and T = 0.3. 
It could be that the case at T = 0.5 is quite sensitive 
to the singularity, and thus it is a somehow intermediate 
case. A more detailed study, either theoretically or by 
considering intermediate or even lower temperatures for 
smaller packing fractions, will be helpful for clarifying 
this issue. On the other hand, the establishment of the 
existence of the two glasses along the line appears to be 
definite by these results. To support this statement, we 
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have plotted in Fig. [l5] both the (partial) static struc- 
ture factor Saa (?) (rescaled by a factor of 2 for having a 
better visualization of the figure) and the non-ergodicity 
factor, respectively at the highest, T — 2.0, and at the 
lowest, T = 0.3, temperatures studied in this work, so to 
compare the most extreme cases of repulsive and attrac- 
tive glasses. Firstly, we note how the oscillations of the 
non-ergodicity factors follow the ones in the structure fac- 
tor quite closely in both cases (for the T = 0.3 case, the 
scale of the figure almost flattens these oscillations, which 
are more evident in Fig. |l4|). Also, the S(q) presents the 
typical features we expected from theoretical calculations 
within MCT and integral equations. Indeed, the repul- 
sive case shows a marked first peak, which is the main 
responsible for the glass transition, while the attractive 
one possesses larger secondary oscillations, that consti- 
tute a signal of the smaller cages already described above. 
Thus, this contributes to establish not only the existence 
of the two glasses, but also, and most importantly, the 
two distinct mechanisms which drive the glassification in 
the two cases, i.e. simply packing and localised attrac- 
tion. 



VI. CONCLUSIONS 

In this paper, using molecular dynamics, we have stud- 
ied the dynamical arrest phenomena of spherical particles 
interacting via a square well potential. The square well 
potential has been studied as one of the simplest canon- 
ical models of solids liquids and gases for many years 
p8| p3[ . Here we have extended the models applicabil- 
ity to the domain of dynamical arrest, and glassy phe- 
nomena. Previous predictions for dynamical arrest from 
mode coupling theory for the square well potential are 
available [|JH|, so direct comparisons are feasible. 

By using a well-adjusted binary mixture, we have been 
able to extent our previous preliminary investigations fl7| l 
much closer to the arrest transition, accessing diffusion 
constants that are 4 orders of magnitudes smaller than 
in the previous calculations. Nevertheless, results on the 
lowest valued iso-diffusivity curve available for the single- 
component system are very close to those for the binary 
mixture, so we may consider the role of the second com- 
ponent to be mainly the prevention of crystallization. 

In that regime where repulsions dominate, we recover 
an ideal-glass-transition with power law scaling of the dif- 
fusion constant. We also find an attractive branch to the 
dynamical arrest where theory has predicted the presence 
of an 'attractive' glass. Where attractions and repulsions 
compete in nearly equal terms we find re-entrance in the 
arrest curves when plotted in units of the microscopic 
temperature dependence of the diffusion constant. How- 
ever, fixed-density diffusion constants, plotted without 
any normalization, exhibit a maximum in the diffusion 
constant as temperature is increased. This maximum 



locates the re-entrant liquid where mobility is anoma- 
lously high. We consider that the essential features of 
re-entrance in this regime for the square-well potential 
now to be confirmed, in agreement with the predictions 
of theory ||]. 

We have studied also the evolution of the density- 
density correlation functions (dynamical structure fac- 
tors) and, independently, the mean-squared distance 
traveled by particles in the vicinity of the reentrant 
regime. As expected, where repulsive interactions dom- 
inate, we find the classical arrest scenario in which 
plateaus develop in both functions as arrest is approached 
These plateaus indicate the development of an ob- 
servable characteristic cage time, and are quite typical 
of prediction by MCT for hard sphere systems. When 
attraction begin to compete on equal terms, in the vicin- 
ity of the re-entrant regime, the theory has predicted the 
existence of an A3 singularity embedded in the arrested 
region. It is therefore difficult to access this singular- 
ity directly by molecular dynamics, but the theory has 
indicated that there are distinctive signatures of this sin- 
gularity in the re-entrant fluid phase, on approach to ar- 
rest. In particular, density correlators from a suitable 
fixed-temperature cut of the phase diagram have an in- 
teresting pattern of behavior in which the logarithmic 
behaviour |13|, due to the embedded A3 point, first be- 
gins to dominate, and then crosses over to the conven- 
tional Ai singular behavior more commonly observed for 
normal MCT arrest. The density correlators in the re- 
entrant regime clearly exhibit this phenomenon, the pat- 
tern of evolutions being essentially in agreement with the 
predictions of theory. 

We may pause here to comment that we do not re- 
gard this complex cross-over behavior as a complication, 
but in fact as a rather delicate, and unusual signature of 
the whole re-entrant phenomenon, and interplay between 
A2 and A3 singularities. That the simulations would re- 
produce this is strong support for the detailed picture 
offered by theory. The behaviour of the mean-squared 
displacements are also quite unusual, and there is as yet 
no theoretical prediction for them in this regime. 

Finally, we are able to extract the non-ergodicity 
factors along the arrest curve, for the lowest iso- 
(normalised)-diffusivity constant curve. Again, in line 
with theoretical predictions we find strong evidence of a 
transition from a repulsive glass to attractive glass behav- 
ior, as indicated by the change in characteristic shapes 
of the non-ergodicity factors. This is the first direct evi- 
dence of the repulsive glass-to-attractive glass transition 
that has been predicted by the theory, representing one 
of the most remarkable phenomena associated with the 
system. 

The theory suffers from strong systematic shifts of all 
the arrest curves in relation to the simulated ones, a 
phenomenon long known from the example of the hard 
sphere. However, in detail the theoretical predictions 
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of re-entrant regime, with associated crossover to loga- 
rithmic singularity, and glass-to-glass transition has been 
confirmed by detailed molecular dynamics calculations. 

From the experimental point of view, there is accumu- 
lating evidence that all the phenomena described here 
are robust, being relatively independent of the details of 
the experimental system used to study them ]l5|,|l6| ,|9| . 
The same is true of the theoretical studies fi|^,[iH7j and 
simulations |s|jL7|]. The square well potential is one of 
the simplest examples one can study, and it is reassuring 
that it exhibits the phenomena. 

Our original prediction that in the dense regime, col- 
loidal particulates system with short-ranged potentials 
could be described using ideas from dynamical arrest and 
glass theory now seems to be strongly supported. Dense 
particle gels are thereby identified as an example of a new 
type of glass, or dynamically arrested phase. The impli- 
cations are broader than the simple example studied, for 
it indicates that it may be possible to interpret many for- 
merly disparate phenomena such as coagulation, precip- 
itation, aggregation, and gellation within the paradigm 
of dynamical arrest or glass theory. This is a fundamen- 
tal soft of perception in the field of dense soft-condensed 
matter which may prove to be very fruitful in coming 
years. 
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observed in Fig. [H](b), because we are here looking at a 



10 



point slightly less close to the glass transition, thus local- 
ization is a bit less effective for this case. However, when 
this happens, it is always inside the attractive well in this 
part of the phase diagram. 
[39] Actually, using the stretched exponential, one can fit the 
correlators not too badly. However, the parameters that 
are found from the fits are completely unreliable. Indeed, 
for the amplitude parameter, which should give an esti- 
mate of the non-erogdicity parameter, we find for every 
q- vector, the value of 1. This is the reason why we do not 
consider these fits to have any validity. 



1 * X i X K 



0.75 

's AA (q*) 

0.5 h 



0.25 



One-component 

X Binary system 



o 
to 




to 



FIG. 1. Comparison between the density-density cor- 
relation function Saa(q*, t)/ Saa{<1*) in the binary and the 
monodisperse case at <j> = 0.57 and T = 0.75, corresponding 
to the most reentrant point found in the mono-disperse case, 
before crystallization intervenes. The wave-vector chosen cor- 
responds for both cases to q" = 2-k/oa- 
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FIG. 2. Normalised diffusion costant D/Dq, with 
Do = ay^T/m in function of packing fraction cj>, along each 
studied isotherms between T = 0.3 and T = 2.0. The normali- 
sation factor takes into account the difference due to different 
inital velocities, and ensures the common low density limit 
(see 0). 




FIG. 3. As in Figj^ but, in this case, the diffusion costant 
has not been normalized, and it is plotted against temperature 
along isochores between <j> — 0.56 and (j> — 0.60. The maxi- 
mum in the bare diffusivity becomes more evident (almost 2 
orders of magnitude) as one moves in the more reentrant part 
of the (0, T)-pfane. 
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FIG. 4. Curves of iso-normalised-diffusivity D/Do in the 
(4>,T) plane. We indicate with a vertical arrow the location 
for the hard sphere case of the packing fraction (<f> ~ 0.582) 
for the lowest iso-D/Do curve, i.e. D/Do — 5e — 6. The inset 
shows the MCT prediction, calculated in ji| , with as input the 
structure factor obtained using Percus-Yevick approximation. 
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FIG. 5. a) Normalised density-density correlation function 
along the isotherm T = 2.0 for different packing fractions. 
The g-vector displayed in this picture is slightly larger then 
the one for the first peak of the static structure factor. The 
horizontal line represents the correspondent non-ergodicity 
parameter, as extrapolated from the fit (see text); 
b) Mean square displacement along the isotherm T = 2.0 for 
the same densities as in figures 0(a) . 
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FIG. 7. a) Same as in Fig.|(a) for T = 0.6. The dashed 
lines represent fits with logarithmic laws. They are displayed 
to show the presence of a logarithmic decay and the mecha- 
nism of its disappearance in the proximity of an A2 singularity 
(see text for further details). They are reported as a guideline 
to the eye, and not to extrapolate any fit parameters, b) As 
in Figj|(b) for T = 0.6. 
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FIG. 8. a) As in Fig.|(a) for T = 0.5. The dashed line is a 
fit of one of the correlators with a stretched exponential, which 
we used to extrapolate the non-ergodicity parameter f q , and 
it is shown to display the goodness of the fit. Parameters of 
the stretched exponential fits are reported in Fig.^. b) As in 
Fig.|(b) for T = 0.5. 
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FIG. 9. a) Exponent f3 q as a function of q obtained 
from the fit of the density-density correlation function 
with the stretched exponential in Eq.|| for temperatures 
T — 0.6, 0.5, 0.3. The values found are always very small, in- 
dicating a very slow relaxation, b) As in a) for the relaxation 
time parameter r q of Eq.^| Interestingly, a peak correspond- 
ing to the q— value of the static structure factor develops, as 
one lowers the temperature. 
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FIG. 10. a) As in Fig.|(a) for T = 0.4. b) As in Fig.|(b) 
for T — 0.4. In the inset the fit of the sub-diffusive and 
subdiffusive regime with a power law is shown (see text for 
details) . 
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FIG. 11. a) As in Fig.|(a) for T = 0.3. In the inset fits 
with streched exponential are shown for different value of q 
at the same packing fraction <f> — 0.58. b) As in Fig.^|(b) 
for T = 0.3. In the inset a fit as in FigJlo|(b) is shown. There 
is no evident subdiffusive regime. 
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FIG. 12. Density-density correlation function along 
iso-D/Do line D = 5 • 10 _6 .The wave- vector chosen corre- 
sponds for all cases to q* = 2n/aA- 
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FIG. 13. Mean square displacement along iso-D/Do line 
D = 5'10" 6 . 
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FIG. 15. Non ergodicity parameter f q and partial static 
FIG. 14^ Non ergodicity parameter /, along iso-D/D line structure factor S g at T = 2.0 and T = 0.3. The different 
D = 5 • 10 6 as obtained by fitting the correlation functions shape of f q re flects the difference between attractive and re- 
(see text for details). pulsive glass (see text for details). 
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